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Abstract 

Mutually unbiased bases (MUBs), which are such that the inner product between two vectors in 
different orthogonal bases is a constant equal to 1/ y/d, with d the dimension of the finite Hilbert 
space, are becoming more and more studied for applications such as quantum tomography and 
cryptography, and in relation to entangled states and to the Heisenberg-Weil group of quantum 
optics. Complete sets of MUBs of cardinality d + 1 have been derived for prime power dimensions 
d = p m using the tools of abstract algebra. Presumably, for non prime dimensions the cardinality 
is much less. 

Here we reinterpret MUBs as quantum phase states, i.e. as eigenvectors of Hermitean phase 
operators generalizing those introduced by Pegg & Barnett in 1989. We relate MUB states to 
additive characters of Galois fields (in odd characteristic p) and to Galois rings (in characteristic 
2). Quantum Fourier transforms of the components in vectors of the bases define a more general 
class of MUBs with multiplicative characters and additive ones altogether. We investigate the 
complementary properties of the above phase operator with respect to the number operator. We 
also study the phase probability distribution and variance for general pure quantum electromagnetic 
states and find them to be related to the Gauss sums, which are sums over all elements of the field 
(or of the ring) of the product of multiplicative and additive characters. 

Finally, we relate the concepts of mutual unbiasedness and maximal entanglement. This allows to 
use well studied algebraic concepts as efficient tools in the study of entanglement and its information 
aspects. 
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I. INTRODUCTION 



In quantum mechanics, orthogonal bases of a Hilbert space 7i q of finite dimension q are 
mutually unbiased if inner products between all possible pairs of vectors of distinct bases 
are all equal to 1/y/q- Eigenvectors of ordinary Pauli spin matrices (i.e. in dimension q = 2) 
provide the best known example. It has been shown that in dimension q = p m which is 



the power of a prime p, the complete sets of mutually unbiasec 
Fourier analysis over a Galois field F q (in odd characteristic p) 



even characteristic 2)|2(. In^O], one can find an exhaustive literature on MUBs . Complete 
sets of MUBs have an intrinsic geometrical interpretation, and were related to discrete phase 



bases (MUBs) result from 
1] or of Galois ring R Am ( m 



spaces 
2- designs 



6], finite projective planes 



Q, 3], convex polytopes and complex projective 
i relation tosvi 



id, There are hints on the relation t o sy mmetric informationally complete 

and to Latin squares 161 ] . 



positive operator measures (SIC-POVMs)[ 

There are strong motivations to embark on detailed studies of MUBs. First, they en- 
ter rigorous treatments of Bohr's principle of complementarity that distinguishes between 
quantum and classical systems at the practical level of measurements. This fundamental 
quantum principle introduces the idea of complementary pairs of observables in the sense 
that precise measurement of one of them implies that possible outcomes of the other (when 
measured) are equally probable. In the nondegenerate case, if an observable O represented 
by a q times q hermitian matrix is measured in a quantum system prepared in the eigenbase 
of its complementary counterpart O c , then the probability to find the system in one of the 
eigenstates of O is just 1/q as corresponding to mutually unbiased inner products. Another 
domain of applications where MUBs have been found to play an important role is the field 
of secure quantum key exchange (quantum cryptography). In the area of quantum state 
tomo grap hy, one should use MUBs for a complete reconstruction of an unknown quantum 
state |l7|. 

In this paper we approach the MUBs theory from the point of view of the theory of 
additive and multiplicative characters in Galois number field theory. The multiplicative 



characters ipk{ n ) — exp 



1 2iirnk \ 
g-1 )■> 



k = 0...q — 2, are well known since they constitute the basis 



for the ordinary discrete Fourier transform. But in order to construct MUBs, the additive 
characters introduced below are the ones which are useful. This construction is implicit in 
some previous papers Q, Q, |^, and is now being fully recognized 13] • 
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An interesting consequence is the following: the discrete Fourier transform in Z q has 
been used by Pegg & Burnett j^J as a definition of phase states \6k), k = 0...q — 1, in TC q . 
The phase states \9k) could be considered as eigenvectors of a properly defined Hermitian 
phase operator Op#. Phase properties and phase fluctuations attached to particular field 
states were extensively described. In particular the classical phase variance ir 2 /3 could be 
recovered. 

We construct here a phase operator Goal having phase MUBs as eigenvectors. In contrast 
to the case of 0p#, we find that the phase fluctuations of Gcai can be expressed in terms of 
Gauss sums over the finite number field F q , and could be in principle smaller than those due 
to Qpb- This points to the fact that the phase MUBs may be of interest for quantum signal 
processing. Character sums and Gauss sums which are useful for optimal bases of m-qudits 
(p odd) are also generalized to optimal bases of m-qubits (jp — 2). 



II. PHASE MUBS IN ODD PRIME CHARACTERISTIC 



A. Mathematical preliminaries 

The key relation between Galois fields F q and MUBs is the theory of characters. This 
has not been recognized before and here we use the standpoint of characters as the most 
general way of considering previous results and also as a better criterium for elaborating on 
future results. 

A Galois field is a finite set structure endowed with two group operations, the addition 
"+" and the multiplication The field F q can be represented as classes of polynomials 



obtained by computin 
the integers modulo p 



? modulo an irreducible polynomial over the ground field F p = Z p , 
2l\ . A Galois field exists if and only if q = p m . We also recall that 
F q [x] is the standard notation for the set of polynomials in x with coefficients in F q . 

A character n(g) over an abelian group G is a (continuous) map from G to the field of 
complex numbers C of unit modulus, i.e. such that \n(g)\ = 1, g e G. 

We start with a map from the extended field F q to the ground field F p which is called 
the trace function 

tr{x) =x + x p + --- + x pm ~ 1 G F p , \/xeF q . (1) 
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Using (|TJ), an additive character over F q is denned as 



K[X I 



exp( 



(2) 



The main property is that it satisfies k(x + y) — k(x)k(u), x, y G F q . 
On the other hand, the multiplicative characters are of the form 



if; k (n)=w£ 1 , k = 0...q-2, n = 0...g-2. 



(3) 



In the present research, the construction of Galois phase MUBs will be related to character 
sums with polynomial arguments f(x) also called Weil sums 



D 



w f = «(/(*))■ 



(4) 



xGF„ 



In particular, ( theorem 5.38 in 2l|), for a polynomial fd{x) G F q [x] of degree d > 1, with 
gcd(d, q) = 1, one gets < {d — l)q 1 / 2 . 

The quantum fluctuations arising from the phase MUBs will be found to be related to 
Gauss sums of the form 

G(V>,«)= £>(aO«(x) , (5) 

where F* = F q — {0}. Using the notation ip for a trivial multiplicative character ip = I, and 
fto for a trivial additive character k = 1 the Gaussian sums (0) satisfy G(ipo,Ko) = q — 1; 
G(ipo, k) = —1; G(?/>, «o) = and |G(/0, k)| = 5 1,/2 for nontrivial characters k and '0. 



B. Galois quantum phase states 

We now introduce a class of quantum phase states as a "Galois" discrete quantum Fourier 
transform of the Galois number kets 

\e^) = ^-J2Mn)<yn)\n), y E F q (6) 

in which the coefficient in the computational base {|0), |1), • • • , \q— 1)} represents the product 
of an arbitrary multiplicative character T/>&(n) by an arbitrary additive character n(yn). 

It is easy to show that previous basic results in this area can be obtained as particular 
cases of (JBJ)- Indeed: 
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Pegg Sc Barnett (1989): For k = k,q and ip = ipk{ri), one recovers the ordinary quantum 
Fourier transform over Z q . It has been shown 20] that the corresponding states 

i**> = 4j£^( n )i n >' (7) 

V q nez q 

are eigenstates of the Hermitian phase operator 

0PB = #fc|0fc)(0fc|> (8) 

with eigenvalues 9k = O + @o an arbitrary initial phase. 

Wootters & Fields (1989): We recover the result of Wootters and Fields in a more general 
form by employing the Euclidean division theorem (see theorem 11.19 in 22|) for the field F q , 
which says that given any two polynomials y and n in F q , there exists a uniquely determined 
pair (a, b) e F q x F g , such that y = an + b, deg(b) < deg(a). Using the decomposition of the 
exponent in (JHJ), we obtain 

|0£) = — ^2 ^ki.n)K{an 2 + bn)\n), a,b <G F q . (9) 

(The result of Wootters & Fields corresponds to the trivial multiplicative character ipo = !)• 
Eq. Q defines a set of q bases (with index a) of q vectors (with index b). Using Weil sums 
(HJ) it is easily shown that, for q odd, so that gcd(2, q) = 1, the bases are orthogonal and 
mutually unbiased to each other and to the computational base 

1 j. ,i \ 2 1 / j l\ $bd if c = a (orthogonality) 

|(^|^)| = |- Y^J;^ c ~ a)n +{d - b)n \ = I k & yj (10) 

1 neFq j -1= if c 7^ a (unbiasedness). 



III. QUANTUM FLUCTUATIONS OF PHASE MUBS IN ODD PRIME CHAR- 
ACTERISTIC 

Following Pegg and Barnett, a good procedure to examine the phase properties of a 
quantized electromagnetic field state is by introducing a phase operator and this was one 
of the reasons that led them to introduce their famous Hermitian phase operator 0p#. In 
Section 6 of their seminal paper they showed "for future reference" how their phase operator 
could be employed to achieve this goal. In this section we proceed along the same lines using 
the phase form of the Wootters-Field MUBs. 
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A. The Galois phase operator 

On the other hand, the phase MUBs as given in (JHJ) are eigenstates of a "Galois" quantum 
phase operator 

0Gai = ]T«>(^|, a,beF g . (11) 

bdF q 

with eigenvalues 9^ = — . We use this fact to perform several calculations of quantum phase 
expectation values and phase variances for these MUBs. 

Using (jHJ) in (fTTjl and the properties of the field theoretical trace the Galois quantum 
phase operator reads 

Ocai = % Yl ~ m)< [a(n2 - m2)] S(n,m)|n)H , (12) 

^ m,n£F q 

where S(n,m) = YbeFq^p ■ ^ n ^ ne diagonal matrix elements, we have the partial 

sums 

S(n,n) = ^l, (13) 

so that (n|BGai|^) = 7r( - l? ~ 1 ' > ■ In the non-diagonal matrix elements, the partial sums can be 
calculated from 

1 — x q qx q n xq 

ox = x{i -|- zx -|- ox - -|- • • • -|- qx' ) = x [- 

b&F q 



Ybx b = x(l + 2x + 3x 2 + --- + qx q - 1 )=x[-^-^--^—} = , (14) 

^— ' L (1 — xY 1 — £ x — 1 



where we introduced x = u P and we made use of the relation x q = 1. Finally, we get 

S(m,n) = - |^y- (15) 

1 - UJp 

B. The Galois phase- number commutator 
Using (fT2|) and the Galois number operator 

N = Ym(ii (is) 

lGF q 

the matrix elements of the phase- number commutator [0Gai, N] are calculated as 

u Ga i(n,m) = ^(n-m)Mn-™K rHn2 ~ m2)] S(n,m). (17) 

The diagonal elements vanish, the corresponding matrix is anti-Hermitian since UGai{n, m) = 
-MQ al (m, n), and the states are pseudo-classical since lim^oo w Gal (n, m) = 0. These prop- 
erties are similar to those of the Pegg & Barnett commutator. 



C. Galois phase properties of a pure quantum electromagnetic state 



For the evaluation of the phase properties of a general pure state of an electromagnetic 
field mode in the Galois number field we proceed similarly to Pegg & Barnett. Thus, we 
consider the pure state of the form 

\f) = X] witn Un = — F ex P(^/3), (18) 

neF q V q 

where (3 is a real parameter, and we sketch the computation of the phase probability distri- 
bution | < 9 b \f > | 2 , the phase expectation value < Ocai >= J2b€F q < ®b\f > | 2 and the 
phase variance < A0Q al >= ^ feeF? (#fe— < ©Gai >) 2 | < @b\f > \ 2 , respectively (the upper 
index a for the base is implicit and we discard it for simplicity). 
The two factors in the expression for the probability distribution 

ipkj—n) exp(m/3)/t(— an 2 — bn)][ V ; fc( m ) exp(— im(3)K(am 2 + bra)], (19) 

n&F q meF q 

have absolute values bounded by the absolute value of generalized Gauss sums G(ip, k) = 
^2x&F q ' l l J (9( x )) K '(f( x ))^ with Fq[ x \- Weil j^j showed that for f(x) of degree d with 

gcd(d, q) = 1 as in (@J), under the constraint that for the multiplicative character ip of order 
s, the polynomial g(x) should not be a sth power in F q [x] and with v distinct roots in the 
algebraic closure of F q , the order of magnitude of the sums is (d + v — \)^fq. For a trivial 
multiplicative character ip , and (3 = 0, the overall bound is | < 9b\f > \ 2 < - and it follows 
that the absolute value of the Galois phase expectation value is bounded from above as 
expected for a common phase operator 

I < e Gal > | < 4 b ^ n - ( 2 °) 

q b€F q 

The exact formula for the phase expectation value reads 

<e Gal >=^ V eP(m,n)S(m,n), (21) 
Q „ 

m,n£F q 

where e /3 (m,n) = ipkijn — n) exp[i(ri — m) f3]x[a(rn 2 — n 2 )] and the sums 5(m, n) were defined 
in (|13|) and (|15|). The set of all the q diagonal terms m = n in < 0Gai > contributes an 
order of magnitude "^qS(n, n) ~ ir. The contribution from off-diagonal terms in (J21)) are 
not easy to evaluate analytically; we were able to show that for them one has \S(m, n)\ = 
|| sin[-tr(n — m)]| -1 . 



The phase variance can be written as 

< A9 2 GaJ >= J2( 6 l ~ 26 » < Gai >)l < 0b\f > I 2 - (22) 

bdF q 

The term < 9c a i > 2 J2beF q I < @b\f > 1 2 does not contribute since it is propor- 
tional to the Weil sum Ylb^F ^p*^™ - ™'* = 0. As a result a cancellation of the quantum 
phase fluctuations may occur in (|22|) from the two extra terms of opposite sign. But 
the calculation are again not easy to perform analytically. For the first term one gets 
2(27r/g 2 ) 2 ^ mneF9 e /3 (m,n)|5(m,n)| 2 . The second term in 022J is -2^,^ < ©Gai > 
I < &b\f > | 2 = — 2 < 6cai > 2 - Partial cancellation occurs in the diagonal terms of (J22J) 
leading to the contribution « — ^ which is still twice (in absolute value) the amount of 
phase fluctuations in the classical regime. A closed form for the estimate of the non-diagonal 
terms is still an open problem. 



IV. PHASE MUBS FOR m-QUBITS 
A. Mathematical preliminaries 

The Weil sums (jlj) which have been proved useful in the construction of MUBs in odd 
characteristic p (and odd dimension q = p m ), are not useful in characteristic p = 2, since in 
this case the degree 2 of the polynomial fd(x) is such that gcd(2, q)=2. 

An elegant method for constructing complete sets of MUBs of m-qubits was found by 
Klappenecker and RottelerQ]. It makes use of objects belonging to the context of quaternary 



codes |24j . the so-called Galois rings i? 4 m; we refer the interested reader to their paper for 
more mathematical details. We present a brief sketch in the following. 

Any element y G R4™, can be uniquely determined in the form y = a + 2b, where a and 
b belong to the so-called Teichmiiller set T m = (0,1, £, ••• ,£ 2m-2 ), where £ is a nonzero 
element of the ring which is a root of the so-called basic primitive polynomial h(x) j^. 
Moreover, one finds that a = y 2 ™ 1 . We can also define the trace to the base ring Z± by the 
map 

171 — 1 

tr(y) = Y,Ay), (23) 

fc=0 

where the summation runs over R^m and the Frobenius automorphism o reads 

cr(a + 2b) = a 2 + 2b 2 . (24) 



In the Galois ring of characteristic 4 the additive characters are 



K(X) = UJ 



tr ( x ) _ ,-tr(i) 



The Weil sums (J3J) are replaced by the exponential sums 
which satisfy 

if y e 2T m , y + 
\V(y)\ = { 2 m if y = 

v2^" otherwise. 
Gauss sums for Galois rings wore constructed [3] 



(25) 



(26) 



(27) 



4m . 



(28) 



where the multiplicative character ip(x) can be made explicit 

Using the notation i(jq for a trivial multiplicative character and k~o for a trivial additive 
character, the Gaussian sums (|25|) satisfy G v (iI)q,Kq) = 4 m ; G y (ip,Ko) = and |G 3/ (^,ft;)| < 
2 m . 



B. Phase states for m-qubits 

The quantum phase states for m-qubits can be found as the "Galois ring" Fourier trans- 
form 

|0<»>) = -L Y, Mn)k(yn)\n), y G R^. (29) 

V neT m 

Using the Teichmiiller decomposition in the character function k one obtains 

|««) = -i^ fc (n)K[(fl + 2%]|n), a,6eT m . (30) 

V neT m 

This defines a set of 2 m bases (with index a) of 2 m vectors (with index b). Using the 
exponential sums (|26jl. it is easy to show that the bases are orthogonal and mutually unbiased 
to each other and to the computational base. The case ijj = ipo was obtained before Q]. 
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C. Phase MUBs for m-qubits: m = 1, 2 and 3 



For the special case of qubits, one uses tr(x) = x in (|30|) so that the three pairs of MUBs 
are given as 

[|0>,|1>]; ^[|0> + |1>,|0>-|1>]; -^[|0>+i|l),|0)-i|l>]. 
For 2-qubits one gets a complete set of 5 bases as follows 

(|0>,|1>,|2>,|3»; 

|[|0) + |1) + |2) + |3), |0) + |1> - |2> - |3), |0) - |1) - |2) + |3), |0> - |1> + |2> - |3>] 
|[|0) - |1) - i\2) - z|3), |0) - |1) + i\2) + i\3), |0) + |1) + i\2) - z|3), |0> + |1) - i\2) + i\3)] 
I[| ) - - i\2) - |3), |0) + i\2) + |3), |0) + i\2) - |3), |0) - i\2) + |3)] 

i[|0) - - |2) - i|3), |0) - i|l> + |2) + i|3), |0) + + |2) - i|3>, |0> + - |2) + z|3>], 

(31) 

and for 3-qubits a complete set of 9 bases 

(|0),|1),|2),|3),|4),|5),|6),|7)); 
i[|0) + |1) + |2) + |3) + |4) + |5) + |6> + |7>, |0) + |1> - |2) + |3) - |4) - |5) - |6> + |7), 
|0) - |1) + |2) - |3) - |4) - |5) + |6> - |7), |0) + |1> - |2) - |3) - |4> + |5) + |6> - |7>, 
|0) - |1) - |2) - |3) + |4> + |5) - |6> + |7>, |0) - |1> - |2) + |3) + |4> - |5) + |6) - |7>, 
|0) - |1> + |2) + |3> - |4) + |5) - |6> - |7>, |0) + |1) + |2) - |3) + |4> - |5) - |6> - |7)], 

(32) 

where only the first two bases have been written down for brevity reasons. 

Quantum phase states of m-qubits (}3T)j) are eigenstates of a "Galois ring" quantum phase 
operator as in (fTTj) . and calculations of the same type as to those performed in Sect, (jlllj) 
can be done, since the tr operator (}2*3*j) fulfills rules similar to the tr operator By analogy 
to the case of qudits in dimension p m , p an odd prime, phase properties for sets of m-qubits 
heavily rely on the Gauss sums (J28|) . The calculations are tedious once again but can in 
principle be achieved in specific cases. 
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V. MUTUAL UNBIASEDNESS AND MAXIMAL ENTANGLEMENT 



Roughly speaking, entangled states in H q cannot be factored into tensorial products of 
states in Hilbert spaces of lower dimension. We show now that there is an intrinsic relation 
between MUBs and maximal entanglement (see below). 

We start with the familiar Bell states 

(\B ,o), |So,i» = ^(|00) + |11), |00) - |11)), (|B 1)0 ), |Bi,i)) = -J=(|01) + |10), |01) - |10» 

where the compact notation 1 00) = |0)©|0), 1 01) = |0) 11),. .., is employed for the tensorial 
products. 

These states are both orthonormal and maximally entangled, i.e., such that 
trace 2 \Bh t k) (&h,k\ = \h, where trace2 means the partial trace over the second qubit j^ . 
One can define more general Bell states using the multiplicative Fourier transform ((7j) 



applied to the tensorial products of two qudits [2 

|B M ) = — J2^ n \n,n + h), (33) 

These states are both orthonormal, (B^By ,k') — 3hh'8kk', and maximally entangled, 
trace 2 \B htk )(B hjk \ = \l q - 

We define here an even more general class of maximally entangled states using the Fourier 
transform @ over F q as follows 

1 q ~ l 

\Kt) = ^Y.< r[{an+b)n] \ n ' n+h ) ■ ( 34 ) 

The h we use here has nothing to do with the polynomial h(x) of Section (JT1J). A list of the 
generalized Bell states of qutrits for the base a = can be found in j3| which is a work 
that relies on a coherent state formulation of entanglement. In general, for q a power of a 
prime, starting from ()34)1 one obtains q 2 bases of q maximally entangled states. Each set of 
the q bases (with h fixed) has the property of mutual unbiasedness. 

Similarly, for sets of maximally entangled m-qubits one uses the Fourier transform over 
Galois rings (}3T)j) so that 

\ B h) = -iE ! * +a,n| l n > n + ^ ( 35 ) 



/Or, 
Z n=0 
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For qubits (m = 1) one gets the following bases of maximally entangled states (in matrix 
form, up to the proportionality factor) 

(|00> + |11>,|00>-|11» (|01> + |10>,|01)-|10» 1 

(36) 

(|00) |00) -z|H» (|01) +i|10), |01>) -i|10)) 

Two bases in one column are mutually unbiased, while vectors in two bases on the same line 
are orthogonal to each other. 

For two-particle sets of quartits, using Eqs. (J3T|) and (JHSJ), one gets 4 sets of (\B^ b ), 
h = 0, 3), see them below, each entailing 4 MUBs (a = 0, 3): 

{(|00) + |11) + |22) + |33) | + + --) I + 1 + - + -)); 

(|oo> — in> — »|22> — *|33> ! + -(+,•)(+•)> l + +(+*)(-0> I ++(-<)(+•*)»;•■■} 



{(|01) + |12) + |23) + |30) | + + --) I + 1 + - + -)); 

(|01)-|12)-i|23)-i|30) | + -(+i)(+i)> | ++ (_i)(+i)));...} 

{(|02) + |i3) + |20) + |3i) | + + --) I + ! + - + -»;•••} 

{(|03) + |10) + |21) + |32) | + + --) I + | +- + -));•••} 

where, for the sake of brevity, we omitted the normalization factor (1/2) and the bases in 
the sets have been labeled by their coefficients unless for the first base. Thus, in the first set 
1 00) + 1 11) + 1 22) + 1 33) = | + + + +). Within each set, the four bases are mutually unbiased, 
as in ([31)1 . while the vectors of the bases from different sets are orthogonal. 

As a conclusion, the two related concepts of mutual unbiasedness and maximal entangle- 
ment derive from the study of lifts of the base field Z p to Galois fields of prime characteristic 
p > 2 (in odd dimension), or of lifts of the base ring Z4 to Galois rings of characteristic 4 (in 
even dimension). One wonders if lifts to more general algebraic structures would play a role 
in the study of non maximal entanglement. We have first in mind the nearfields that are 
used for deriving efficient classical codes and which have a strong underlying geometry [3^. 



VI. CONCLUSION 



In this research, we approached the MUBs fundamental topic from the point of view of 
the additive and multiplicative characters over finite fields in number theory. We consider 
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that this framework is the most general including previous results in the literature as par- 
ticular cases. Since MUBs are essentially generalized discrete Fourier transforms over finite 
number field kets, we formulated a quantum phase interpretation and illustrated several 
calculations of the phase properties of pure quantum states of the electromagnetic field in 
this finite number field mathematical context. Various types of Gauss sums get involved 
in this type of calculations of the MUBs phase properties of a pure quantum state and 
the generalization to the mixed states, although straightforward through the usage of the 
density matrix formalism, could lead to even more complicated calculations involving such 
sums. We hope to evaluate them in future works. We also mentioned in the last section a 
possible application to phase MUBs states of Bell type. This could lead to finite number 
field measures of the degree of entanglement. 

Note: The authors acknowledge Dr. Igor Shparlinski for suggesting important corrections 
on this paper. On June 20/2005 he sent us an e- message pointing to some mathematical 
inconsistencies in subsection III.C (i.e., subsection 3.3 in the epjd version). Parameter (3 
therein should be itself an element in F q to perform the calculations. 
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